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Abstract 

Image segmentation is an important step in texture 
recognition of Satellite Synthetic Aperture Radar (SAR) 
images. In this paper, an efficient method for texture 
recognition of SAR images is proposed based on 
extraordinary result of wavelet transform on texture feature 
extraction and the benefits of genetic algorithm (GA) as a 
classifier. First, texture feature is extracted by wavelet 
transform. Second, a feature vector composed of wavelet 
energy feature, kurtosis value of wavelet energy feature and 
gray values of eight-neighborhood of SAR image is formed. 
Then, segmentation of different textures is applied by using 
feature vector and level set function. Finally, the values of 
feature vector are used for training the classifier. Two 
different classes of datasets are used to show the 
experimental results. One is a simulated SAR image and the 
other is a real satellite SAR image. The results of experiments 
show the powerful ability of this algorithm to recognize 
different textures in Satellite SAR images. 
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Introduction 

Synthetic aperture radar (SAR) imaging has many 
applications for environmental monitoring, recognition, 
targeting, surveillance, earth-resource mapping, and 
military systems. SAR imaging has capabilities to 
provide broad-area imaging at high resolutions and to 
acquire images in inclement weather or during night as 
well as day. These applications are driven by the 
military's need for all-weather, day-and-night imaging 
sensors. Because of the ability to use in all weather 
conditions, day and night xerography time and the 
high spatial resolution, SAR imaging systems are 
known as the most popular remote sensing technique 
widely used in the past decades [1]. In [35], a portable 
chip based on SAR imaging system is designed which 



can be installed on SAR satellite and earth stations 
(terrain receivers) for transform the microwave signal 
to image as shown in Fig. 1. This Fig. describes how a 
Satellite SAR imaging system works. 



SAR Ssitf Uirr 




Eniiti Statin 

FIG. 1 THE CIRCUIT CHIP BASED ON SAR SATELLITE IMAGING 
SYSTEM PROPOSED IN [35] 

In remote sensing applications, SAR imaging is a very 
powerful and attractive tool due to its high spatial 
resolution [2]. However, automatic interpretation of 
SAR images is extremely difficult because of the 
speckle noise [3]. The speckle noise is a perfectly 
expanded noise which SAR images are commonly 
modeled as influenced by it. SAR speckle can be 
described as multiplicative noise, with standard 
deviation equal to pixel reflectivity value [4]. The 
speckle noise also affects the probability density 
function (PDF) of the pixel intensities in SAR images. 
This phenomenon can be expressed by the nonlinear 
intensity inhomogeneity in SAR images [2]. This 
phenomenon causes that traditional image processing 
algorithms are not useful to SAR image processing. As 
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a result of this phenomenon, accurate segmentation for 
SAR images is still a crucial issue. So, we should search 
for methods to be independent from intensity in SAR 
image interpretation. 

Segmentation of SAR images is a fundamental process 
in SAR image interpretation. Due to the existence of 
correlated speckles, the segmentation of SAR images is 
commonly expressed as a very complex work [5]. The 
purpose of SAR image segmentation is to partition an 
image into areas of different characteristics [6]. It is a 
crucial point of SAR images interpretation in a way 
that the SAR image is decomposed in a mosaic of 
uniform areas. Recently, a number of methods have 
been proposed for SAR image segmentation including 
global segmentation techniques [7], [8], fuzzy 
clustering methods [9], spectral clustering algorithms 
[6], [10]-[13], and active contour model and level set [2], 
[5]. 

A global segmentation approach adapted to Gamma 
pdf of SAR images is used in [7] for SAR image 
segmentation. This parametric scheme is based on a 
polygonal grid model with a suppositional unknown 
number of regions. In this approach, the number of 
regions of the partition is estimated by minimizing the 
stochastic complexity. However, when the gray values 
of SAR images are not correctly described by Gamma 
pdf, like in textured speckle images, this technique will 
be failed. Furthermore, since a parametric noise model 
of the segmentation technique is considered, the 
parameters of regions need to be tuned, especially 
when the data differ from Gamma pdf such as in 
textured areas. So, this parametric technique only 
describes a unique class of textured data which can 
lead to analogous limitations to those obtained with 
the Gamma pdf. 

In [8], the same Gamma distributed-based method as 
that in [7] was applied with different noise models. 
Compared with segmentation algorithm proposed in 
[7], this approach has some obvious advantages. It is 
based on a polygonal grid which can have an arbitrary 
topology and whose number of regions and regularity 
of its boundaries are obtained by minimizing the 
stochastic complexity of a quantified version on Q 
levels of the image [8]. Unlike in [7], this approach 
leads to a criterion without parameters that need to be 
tuned by the user. However, the appropriate value of 
Q, which minimizes the number of misclassified pixels 



(NMP), could not be estimated automatically in that 
approach. Moreover, if a texture area of an SAR image 
has other characteristics of parametric noise models 
such as Gaussian or Poisson, this technique may fail 
due to the nonlinear intensity inhomogeneity 
phenomenon mentioned above. 

Recently, different forms of spectral clustering (SC) 
algorithms [6], [10]-[13] were proposed. Compared 
with traditional clustering algorithms, SC has some 
obvious advantages. It can recognize the clusters of 
unusual shapes and obtain the globally optimal 
solutions in a relaxed continuous domain by eigen 
decomposition [6]. However, the complexity of 
calculations exists in these methods and they are 
computationally expensive because of using an affinity 
matrix established by the likeness of each pair of pixels. 
So, the method needs to compute the eigenvectors of 
the affinity matrix. Furthermore, SC requires assigning 
a parameter, namely, the scaling parameter a in the 
Gaussian radial basis function (RBF). Appropriate 
assigning of a is determinant to obtain good 
segmentation results in SC. Unfortunately, it is difficult 
to select the suitable o value, and it is always set 
manually. The improper value of o can deteriorate the 
capabilities because SC is highly sensitive to a, and 
different values of a may cause to extremely different 
results [6]. 

In [14] and [15], the techniques established upon the 
active contour model and level set are proposed for 
segmentation of SAR images. These techniques are 
typically based on optimizing an energy function 
depending on the whole image and they can overcome 
the disadvantages of SC models. However, because the 
corresponding energy function has a local minimum, 
these methods may result in an unexpected state. So, in 
some cases, the level set function in these techniques 
cannot converge to a confident interval. 

However, recent SAR image segmentation models [5], 
[9] tend to rely on intensity homogeneity in each 
region as well as no intention for segmenting the 
special objects in SAR images. On the other hand, 
nonlinearity of intensity inhomogeneity often occurs in 
SAR images as discussed above. These two methods 
cannot address the nonlinear intensity inhomogeneity 
phenomenon. Moreover, these methods need an initial 
curve to be created by user. 

In this paper, we propose a new energy function 
named kurtosis wavelet energy (KWE) in order to 
solve these mentioned problems. This energy function 
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is an extension of that in [2] . KWE energy function will 
be concluded by pursuing the higher order statistics 
(HOS) and fourth order normalized cumulant, named 
kurtosis, implemented by wavelet and level sets. This 
energy function is used as a texture feature in this 
paper for SAR image segmentation application and has 
the following advantages: 1) it extracts more statistical 
characteristics of a specific region from an SAR image 
because the KWE energy formulation derived in this 
paper is one_order greater of cumulant than the energy 
function proposed in [2], and 2) it does not need to 
have an initial contour, which should be initialized by 
a user. So, our algorithm can automatically perform the 
segmentation procedure of SAR images. 

A new region-based active contour model for 
automatic segmentation of SAR images is proposed in 
this approach. We first define the kurtosis wavelet 
energy (KWE) function where it locally measures the 
region textures along the contour and texture features 
on two sides of the contour. So the texture of each 
region, can be characterized by using the energy of 
wavelet coefficients of each region in each sub-band of 
the decomposition process. The obtained texture is 
then formed into level set functions with a fourth order 
normalized cumulant term. Then, the regions and their 
interfaces can be modeled with these level set functions 
and a functional formulation can be calculated on these 
level sets in terms of fourth order normalized 
cumulant. Finally, the optimal automatic segmentation 
based method will be defined on textures by 
minimizing of this functional formulation. 

Recently, Gabor filters are used to extract texture 
features for segmentation [16]-[18], but they have a 
major disadvantage, so that they produce a lot of 
redundancy and thus lots of feature channels. To 
overcome this problem, the structure tensor is 
employed to discriminate textures by many authors 
[19]-[23]. But the structure tensor misses so many other 
useful texture features that it is difficult to acquire 
perfect results [24]. 

In recent and thorough studies, Randen and Husoy 
[25], Haralik [26] and Haralik et al. [27] have reviewed 
and compared different approaches to texture feature 
extraction. Randen' s images, which are fairly hard to 
classify even by eye, are becoming a benchmark for 
assessing different segmentation algorithms [28]-[31]. 
However, these techniques are not suitable for texture 
definition and discrimination of SAR images. The 
segmentation or classification methods by Randen et al. 
are not specifically designed for SAR images textures. 



In other words, the existence data as a ground truth for 
SAR images are different from Randen' s data as a 
benchmark for optical texture algorithms. 

Note that our method, termed as kurtosis wavelet 
energy model, is also related to SWE and LBF models 
which were presented in [2] and [32] respectively. In 
[32], the local intensity average values were derived as 
the minimizers of the energy function in a variational 
formulation. Whereas in this paper, a local texture of 
each region of a SAR image is derived as a minimizer 
of the proposed kurtosis wavelet energy function. 

The remainder of this paper is organized as follows. In 
Section 2, we introduce kurtosis as an effective feature 
for SAR image segmentation. In section 3, derivation of 
classification formulation based on genetic algorithm is 
presented. The implementation and experimental 
results of our proposed algorithm on both simulated 
and real satellite SAR images are given in Section 4. 
Finally, Conclusions are drawn in Section 5. 

Kurtosis Wavelet Energy Formulation 
as a Texture Feature 

In probability theory and statistics, the moment 
generating function of a random variable X is defined 
by 

M x (a)) = E x { e j rox } (1) 

The cumulant generating function of the same variable 
X is given by the logarithm of the moment generating 
function, C x (co)as 

C x (g>) = ln(M x (G>)) = ln[E x {e jo)X }] (2) 

After making a Taylor series expansion in the co 
variable about the cumulant generating 
function, C x (co) which is given by equation (2), can 
obtain as below 

C x (cv) = ]n[E x {e j6>X }] = Yk n ^- (3) 

n=\ 

The coefficient of the Taylor series term co n , multiplied 
by (-j) n , is called the nth-order cumulant of the 
random variable X and it is defined by 

Cx n) -(-j) n 7^rC x (co)| m=0 for n = 1,2,3,- (4) 

The first order and the second order cumulants are the 
mean and the covariance of random variable X 
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respectively. It can be found by computing the first 
and the second order cumulants, C (1) x (co) andc (2) x (co), 
from equation (4) and rewriting them in terms of the 
moments. The first order moment is the mean and the 
second order moment is the covariance (or center of 
mass). These relations are defined by 

H = c«=m«=E{X} (5) 

and 

cov. = c< 2) = m< 2 > - (m^) 2 = E{X 2 } - (E{X}) 2 (6) 
respectively. 

It is reported [2] that whatever the order of cumulant 
as a feature for a SAR image increases, this feature will 
give the more statistical characteristics of a specific 
region from a SAR image. However, implementation of 
higher order cumulants will impose more calculations 
and time consuming. So, we should perform a trade off 
between the higher order and the complexity of 
calculations. In this paper, we propose to use the 4 th- 
order normalized cumulant named kurtosis as a 
texture feature for segmentation of SAR images. We 
show that this selection gives the more statistical 
information of each region also it has logical 
computation and efficient implementation. 

The expression yielding the fourth order cumulant of a 
random variable X from its moments can be derived by 
the relations that exist between moments and 
cumulants. The fourth order cumulant formulation of 
the interested variables is the following 

c^=E{(X-u0 4 } 

= -3[m^ 2) ] 2 -4m^.m^ ( 7 ) 
+ 12[m^] 2 .m^-6[m^] 4 

In the above 4 th order cumulant formulation, we can 
use it as a statistical feature for texture discrimination 
of different regions in SAR images. However, the 
formulation of 4 th order cumulant derived in equation 
(7), is not inherently independent from the size of its 
random variable. On the other hand, in the SAR image 
segmentation approach, we should use a feature which 
it is invariant to image size. In particular, it is shown 
[33] that if a is an arbitrary scalar, then the cumulants 
of the random variables X and aX are related as 

Y = aX o 4 m) =a m 4 m) (8) 

In order to make the 4 th order cumulant formulation 
independent from its random variable, we provide a 



normalization of it. For this purpose, we consider the 
normalized cumulant of order (m,n) associated with 
random variable X as defined by 



r x I 

In the above equation, Cx m ' n) is considered as the m th- 
order cumulant normalized with n th-order cumulant. 
Such formulation can solve the image size variant 
problem and will act as image size invariant. Using 
this definition of cumulant for texture discrimination, 
it directly acts according to that this feature 
(normalized cumulant) is scalar invariant. This means 
that C^' n) = C^ m n) for any nonzero scalar, a. 

Typically, the order parameters m and n are selected as 
integer numbers in a way that m>n. Note that in this 
equation, the denominator should not be zero. So, 

( 2 ) 

we select n=2 to provide a logical choice, since c x is 
covariance and it is always nonzero for any random 
variable. If we select m=4 and n=2 and rewrite the 
normalized cumulant formulation derived in (9), we 
achieve the kurtosis formulation which we propose in 
this paper as [1] 

E{(X-|i) 4 } Cl 4) 
kurtosis(X) = / = (10) 
c (var) 

where, \i is the mean of random variable X, var is the 
variance a 2 of random variable X and it is a nonzero 
parameter (var > 0), and E{t} represents the expected 
value of the quantity t. 

In equation (10), Cx' 2) is named kurtosis which is the 
fourth order cumulant normalized with covariance. 
The kurtosis cumulant indicates the poignancy of the 
PDF histogram of a SAR image. In other words, this 
feature is the slope decrease curve of the PDF 
histogram of a SAR image. With a view to statistical 
calculations, kurtosis is a measure of how outlier-prone 
a distribution is. The kurtosis of the normal 
distribution is 3. Distributions that are more outlier- 
prone than the normal distribution have kurtosis 
greater than 3, and distributions that are less outlier- 
prone have kurtosis less than 3. 

We use the kurtosis as a parameter of the Energy 
distribution of the wavelet coefficients in each sub- 
band of a SAR image. 
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Derivation of classification Formulation 
based on genetic algorithm 

As we know, moments with higher orders can 
characterize textures in images. Since cumulants are 
able to give more statistical information, we can use 
high order cumulants to describe a texture of each 
region in SAR images. Thus, a texture discrimination of 
each region in SAR images can be done by the kurtosis 
explained in previous section. 

We consider that a texture of each region in a SAR 

image is characterized by the kurtosis energy of its 
wavelet coefficients. Let T (x,y) be a texture function 
of each region of a SAR image as below [1] 



X Z! W 4> ( j° ' m ' n) * io.m,n O' y) 



(11) 



T (x,y) = - 



VMN 

-i 

Z EZZ w i ( j' m ' n) ^Unfey) 

| i=H,V,Dj=-j m n 

Vmn 

where \\) is the wavelet function, (j) is the scaling 
function and j is the order of the decomposition. So, 
the energy of the wavelet coefficients in each sub-band 
can be obtained by the following sequence 

(| ( j , m, n)| 2 , | w; ( j, m, n)| 2 , m, n e Z,- j < j < 1) (12) 

Aujol et al. [34] developed this model and showed 
experimentally that the distribution of the energy of 
the wavelet packet coefficients in a sub-band of any 
image follows a picky Gaussian form as below [1] 



Px2(y)= ^ exp( " ( 4 ),1) 



|y>0 



(13) 



where the parameters r\ and C were defined to classify 
the textured images. Such model can not classify 
regions of SAR images. In equation (13), rj and C were 
computed from the first and second order moments of 
the energy distribution in each sub-band. Such 
definition is not good for our SAR image segmentation 
problem because of the instability that exists in the 
curve of moments. In our segmentation approach, rj 
and C are the segmentation parameters and we are now 
in position to calculate these segmentation parameters 
in relation to the higher order moments (third and 
fourth order moments) of the energy distribution of the 
wavelet coefficients in each sub-band. 

Suppose that for a given point p □ r, where r is a 
desired region of a SAR image, the kurtosis wavelet 



Energy function F^ WB is followed by the distribution 

form of the kurtosis energy of the wavelet coefficients 
in a sub-band of each region. From equations (7) and 
(10), we have [1] 



C (4) 



Kurtosis = - 



/p(2) \4/2 

_m£>-3[m£>] 2 -4mg>jn£> 

/p(2) \4/2 

12[m^] 2 .m^-6[m«] 4 

+ (p(2) \4/2 

The Gamma function T(t) on R + is defined as 

00 

r(t) = jVv 4 du (is) 

o 

We assume that p ^ 2 (y) defined in equation (13), is a 
probability density function of a sub-band of each 
region. Also, I W I is used instead of the random 

2 2 

variable X . I W I is a symbol to show the energy of 
wavelet coefficients in each sub-band. So, we have 



(14) 



|P,w^ W ») dW »= N 



(16) 



where N is the total number of pixels of the given 
detail SAR image. The value of parameter A in 
equation (13), can be obtained by using the equations 
(13), (15) and (16), and making the change of 

variable^^L^ _ u as below 



A = 



Nr) 

§r(-) 



(17) 



Now, we can calculate the formulations of the first, 
second, third and fourth order moments of wavelet 
coefficients energy in terms of T(t) from equations 

+00 

= E f w * |2 J = I 1 W * |2 w(' W * l)d c W * l) 



— 00 
+ 00 



= J (W^ 2 Aexp(-(-^)^d(Wp (18) 



r| r| 

nc ; 2=E {IW (p l 4 } = — T(-), 



" iwcpr 



(19) 
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and 

<,=E{IW/} = ^V) 

respectively. 

Dividing (21) by (20), we obtain 



(4) 


7 


r(-) 


m iw/ 




(3) 


9 


m iw/ 


r(-) 







(22) 



Also, 7] can be obtained by rewriting equation (14) as 
below 

r (4) 

r, = F (Kurtosis) = F ( (2) 9 3/2 ) (23) 
(m iw/) 

where F(x) is the function obtained by rewriting the 
kurtosis formulation in equation (14) in terms of the 
first, second, third and fourth order moments 
equations derived in equations (18)-(21) as follow [1] 

19 15 13 7 

r 3 (-)r(-) - 3.N.r 2 (-).r 2 (-) - 4.N.r 2 (-).r(-).r(-) 
v/ , LJ ii ii ii ii ii 

W " 1 5 3 

N.(r(-).r(-)-Nr 2 (-)) 2 
ii ii ii 

13 5 3 

i2.N 2 .r(-)r 2 (-).r(-) - 6.N 3 .r 4 (-) 
11 11 11 !_ 

N.(r(-).r(-)-Nr 2 (-)) 2 
ii ii ii 

(24) 

The graph of F ~ 7 (Kurtosis) is shown and compared 
with the other proposed methods in Fig. 2 [1]. 




[2n d Mome nt]o r[Skew n essjo r[K u rtos is] 



FIG. 2 THE CURVE OF F _1 (KURTOSIS), F _1 (SKEWNESS) AND 
F" 1 (2nd MOMENT) [1] 



As depicted in Fig. 2, the curve of F ' (Kurtosis) is more 
stable than the F " 7 (skewness) and F " 7 (2 nd Moment) 
which were proposed in [2] and [34] respectively [35]. 
So, kurtosis is a good texture feature for satellite SAR 
image recognition. Now, we want to use the kurtosis 
values of the energy of wavelet coefficients of each 
texture in SAR images as a unique property to form a 
feature vector in the next section. In order to do that, 
we should study the active contour models. 

Active contour models were first formulated by 
Mumford and Shah for the image segmentation 
problem [36]. Their function is practically difficult to 
be minimized because of the unknown contour and the 
no convexity of the function. Afterwards, Vese and 
Chan [37] optimized the active contour of the 
Mumford-Shah problem for a special case where the 
image is considered as a piecewise constant function. 
Their model does not contain any local intensity 
information. Then, Chunming Li et ah [32] proposed 
the LBF model to optimize the Vese-Cahn model. 
However, their model does not contain any statistical 
information, which is crucial for segmentation of SAR 
images with nonlinear intensity inhomogeneity. In this 
section, we introduce the level set based kurtosis 
model instead of the LBF model. 

In [35], we defined a segmentation problem based on 
any textural region in SAR images. In that work, we 
defined the region Re k = {p e r I p belongs to the region kj 

where r was an open subset of R 2 , K was the number 
of segmented regions (number of r k ), k was a 
parameter that shows the region index, p was a point 
of region, R kl was the interface between region r k and 
region r x and the image was a function considered as I : 
r — > R . Also, we denoted that for all k= 2,..., K, Re k is 
an open set r k given by a Lipschitz function & k ;r— ► R 
so that ® k (p) > if pixel p is inside the region, 
O k (p) = if pixel p lies in the border of the region, 
and otherwise O k (p) < 0, where & k is the signed 
distance function to the region R k . We can also 
determine r k using the sign of & k and the Heaviside 
distribution function H [1]. 

We define the following functional 

Trwe-ga £ KWE-GA 

(flO + HPflD) (25) 
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where \i is a positive constant, s^ WE _ GA is the kurtosis 
wavelet energy based on GA, and p(O) is the deviation 
of the level set function & from a signed distance 
function as below 



p(d>) = jl(|V<D(p)|-l) 2 dp 



(26) 



We propose the kurtosis of the wavelet coefficients 
energy based on GA to be used as contour energy 

s p 



We consider that for a given point per, the kurtosis 
wavelet energy 8^ WE GA is followed by the distribution 
form of the kurtosis energy of the wavelet coefficients 
in a sub-band of each region. 

To find the boundary of regions in SAR images, we 

should minimize the integral of SKWE GA° ver a ^ center 
points p in the all regions of a SAR image. So, we 
define the energy functional s(C,k 1 ,k 2 J as bellow: 

s(C, k x (p), k 2 (p)) = J s p KWE _ GA (C, k x (p), k 2 (p))dp (27) 

r 

The values of the functions k { (p) and k 2 (p) that 
minimize the KWE energy e^E.^ (C, k x , k 2 )are 
functions of the center point p of the kernel function 
K(p-q). Thus, we claim that KWE energy SKWE GA nas tne 
localization property. 

If we represent the contour C cz r with the level set of a 
Lipschitz function &: r — > R, we get 

(0,k 1? k 2 ) = JsP WE _ GA (0,k 1 (p),k 2 (p))dp 

r 

= Yi{[{k |w ^ |2 (p - q) I I(q) - k,(p) I 2 H(0(p))dq]dp 

+ Y2{[{k |w ^,2 (p - q) I I(q) - k 2 (p) I 2 (1 - H(<D(p)))dq]dp 

(28) 

where H is the Heaviside function, and y l and y 2 are 
two constants. 

Now, we define the new entire kurtosis wavelet energy 
based on GA function F KWE _ GA as 

Frwe-ga 

(0,k„k 2 ) 

— KWE-GA 

(O.k^k^ + ^O) (29) 

where \i is a nonnegative constant, and p(O) is given 
by equation (26). 



Keeping k x (p) and k 2 (p) fixed, and minimizing the 
F KWE _ GA (0,k 19 k 2 ), the active contour model proposed 
in this paper, is derived by equation (30) 

— = - 8 (0)(y 1 e 1 - y 2 e 2 ) + v5 (O)div(- —) 
ot VO 



9 VO 
+ Li(V 2 0-div(—-)) 
VO 



(30) 



where the functions e x and e 2 are computed by 
e 1 (p)=fK |w|2 (q-p)II(p)-k 1 (q)l 2 dq (31) 



and 



e 2 (p) = J K |W/ (q - p) I I(p) - k 2 (q) I 2 dq (32 ) 



respectively, and k x (q) and k 2 (q) are given by 
K |w|2 (q)*[H p (<D(q))I(q)] 



and 



k 2 (q) = ^T 



K |w/ (q)*H p (3)(q)) 
K |w|2 (q)*[(l-H p (0(q)))I(q)] 



K |w|2 (q)*[l-H p (<D(q))] 



(33) 



(34) 



Respectively [1]. 



Implementation and Test Results 

To verify the satellite SAR image intensity model and 
the recognition method we proposed, the results of 
different algorithms on simulated and real SAR images 
will be presented in this section. 

Segmentation of Simulated SAR Image 

In this subsection, we first show one typical example of 
our large number of experiments with synthetic 
images and define the measures we adopted for 
performance analysis. We first simulate a 1-look case 
by the synthetic image of Fig. 3. It consists of four 
regions, as visual inspection can quickly indicate. It has 
been obtained by multiplying a perfect reflectivity 
synthetic image by a white exponential speckle noise 
[15]. The corresponding 1-look noisy image, as shown 
in Fig. 3(a), is generated by multiplying a realization of 
speckle. The ground truth image, as shown in Fig. 3(b), 
is used to calculate the error rates of the segmentations 
obtained by different algorithms. 
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(i) (j) 



FIG. 3 (A) ONE-LOOK SIMULATED SAR IMAGE (222 X 222), (B) GROUND TRUTH, (C) INITIAL CURVE, (D) INTERMEDIATE CURVES 
PRODUCED BY KWE-GA METHOD PROPOSED IN THIS PAPER, (E) FINAL POSITIONS, (F) COMPUTED CLASSIFICATION KWE-GA 
METHOD PROPOSED IN THIS PAPER, (G) FINAL CONTOUR EXTRACTED BY THE METHOD PROPOSED IN [15], (H) SEGMENTATION 
OBTAINED BY THE METHOD PROPOSED IN [15], (I) THE DIFFERENCE BETWEEN THE RESULT OF KWE-GA AND GROUND TRUTH 
(ERROR RATE: 0.217%; THE NUMBER OF MISCLASSIFIED PIXELS: 107), AND (J) THE DIFFERENCE BETWEEN THE RESULT OF THE 
METHOD PROPOSED IN [15] AND GROUND TRUTH (ERROR RATE: 0.942%; THE NUMBER OF MISCLASSIFIED PIXELS: 464) 
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Fig. 3(c) shows the initial position of the evolving 
curve, Fig. 3(d) shows an intermediate position of the 
curves, Fig. 3(e) shows final position, and, finally, Fig. 
3(f) displays the segmented and classified regions 
represented by their KWE value at convergence and 
GA to classify. The segmentation results by the method 
proposed in [15] are also shown in Fig. 3 and some 
comparisons between these two methods are reported. 
Fig. 3(g) shows the final contour convergence to 
segment the four regions with the method proposed in 
[15], and Fig. 3(h) shows the segmentation result 
obtained by this method. In comparison between Fig. 
3(f) and 3(h), it can be visually found that KWE-GA 
method is better than the method proposed in [15]. 
Since we have the ground truth, we can also give 
objective evaluation in addition to visual inspection. 
The misclassified pixels and the error rate are 
computed by matching the classification result with 
the ground truth to evaluate the method. Fig. 3(i) 
shows the difference between Fig. 3(f) and 3(b) which 
is matching the KWE-GA result with the ground truth. 
Misclassified pixels are the pixels with non zero values 
in Fig. 3(i). The number of misclassified pixels can be 
estimated by counting the pixels with non zero values 
in Fig. 3(i). The number of misclassified pixels is 107. 
Since total number of pixels in Fig. 3(i) is 
222 x 222=49284, we can compute the error from 
107/49284=0.00217. So, the error rate for the KWE-GA 
method is 0.217%. Finally, Fig. 3(j) shows the difference 
between Fig. 3(h) and 3(b) which is matching the 
method proposed in [15] with the ground truth. The 
number of misclassified pixels and the error rate will 
be computed equal to 464 and 0.942% respectively 
from this Fig. From these reports, we found that the 
overall error rate is reduced about 0.725% by using 
KWE-GA method proposed in this paper. Therefore, 
the KWE-GA algorithm has the best classification 
result according to the error rate. 

Segmentation of Real SAR Images 

In this section, to demonstrate the applicability of the 
KWE-GA algorithm to real SAR image classification, 
experiments on a real SAR image are performed. To 
elucidate relative advantages of the KWE-GA 
algorithm with respect to KWE single method 
proposed in [35], the results of different algorithms on 
a genuine monolook satellite SAR image are presented. 
Not that when we deal with real SAR images 
classification, the ground truth corresponding to the 
SAR images being classified is absent generally. In this 



case, the evaluation of the classification result is based 
on visual inspection of the classified images. 

The SAR image, as shown in Fig. 4(a), is a real single 
band SAR image acquired in August 1989 over the 
agricultural region of Feltwell (U.K.) by a fully 
polarimetric PLC-band NASA/JPL airborne sensor. 
More precisely, all three polarizations (HH, HV, and 
VV) acquired at band C, the HV and VV polarizations 
acquired at band L, and the HH polarization acquired 
at band P were used. The remaining channels (i.e., L- 
HH, P-HV, and P-VV) were discarded, because their 
histograms exhibit strong irregularities (for further 
details on this data set, we refer the reader to [38] -[40]). 
Hereafter, the adopted Feltwell bands will be 
synthetically denoted by // Feltwell-CHH/ , "Feltwell- 
CHV," // Feltwell-PHH. ,, 

To calculate the functions ki(q) and ki(q) given by 
equations (33) and (34) respectively, we should 
compute the kernel function K |W |2 (p) . For this reason, 

we first compute the value of m 2 , m ^ 2 , m ^ 2 and 

IW^I '*V(p' '*^(p' 

m ^ |2 accordingly. To compute the first, second, third 

and fourth order moments of the energy distribution of 
the wavelet coefficients, we use wavelet 
decompositions of the maximum level L. Then, we 
calculate the kurtosis of these coefficients. We have 
tested several kinds of mother wavelets and several 
wavelet coefficients , , and . 

indicates that wavelet coefficients are measured 
along columns (horizontal direction), along rows 
(vertical direction), and along diagonals. , 

and are detail coefficients and contains 

approximation coefficients. We have selected to use the 
Haar mother wavelet as well. 




(a) 
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(b) 



FIG. 4 (A) "FELTWELL-LHV" SAR IMAGE OF U.K. (309 X 433), (B) 
REPRESENTATION OF THE APPROXIMATION WAVELET 
COEFFICIENTS ENERGY | W(p A I 2 FROM LEVEL 1 TO 8 OF THE 

"FELTWELL-LHV" SAR IMAGE 

After calculating the kurtosis of wavelet coefficients 
energy, the value of r\ can be computed from the curve 
of F " 7 (kurtosis) represented in Fig. 2, and then the value 
of C can be obtained from equation (22). 

In Fig. 4(b), we show the approximation wavelet 
coefficients energy I I 2 based on genetic algorithm 

from level 1 to 8 of the "Feltwell-LHV" SAR image. 

The value of the kurtosis of the approximation wavelet 
coefficients energy based on genetic algorithm 
represented in Fig. 4(b) is 1.087 x 10 11 . It is obtained 
by equation (14). Now, we can get t]=0.2837 from Fig. 5. 
This Fig. is the extended curve of Fig. 2 in the range 
[0, 7 x 10 n ] of only the kurtosis-GA axis in detail. 

Also, we calculated C = 0.2932 and 

A= 3.7896 x 10 4 from equations (22) and (17), 

respectively. Note that the kurtosis parameters, r\, C 
and A are different for each SAR image. 
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FIG. 5 THE CURVE OF H IN TERMS OF KURTOSIS OF FIG. 4(B) 

To implement our KWE-GA texture recognition 
algorithm, the parameters in the experimentation are 
supposed as follows: y 1=1.0, yi=1.0, v = 0.004 x255 2 , 
c =2 (constant value of step function used as initial 



contour), time step x-0.1, ju = 1.0, and center point 
p=1.0 . 

On SAR image as shown in Fig. 4(a) which is shown in 
Fig. 6(a) again, three algorithms are used for 
classifying analysis: 1) SWE model, 2) KWE active 
contour model, and 3) the proposed KWE-GA. 

The segmentation obtained by the SWE model [2] is 
shown in Fig. 6(b). One can see that the water area 
(middle top) is incorrectly segmented. Furthermore, 
the boundary between the vegetation and the water is 
not correctly defined (for better evaluation see Fig. 
6(e)). The segmentation obtained by KWE, as shown in 
Fig. 6(c), improves the uniformity in the water region. 
However, there is still mis-segmentation in the water 
region [see the magnified of a selected area in Fig. 6(f)]. 
KWE-GA gets the best classification result, as shown in 
Fig. 6(d). The classification obtained by KWE-GA, 
improves the uniformity in the water region, and a 
local region of the crop in the water [see zoomed image 
in Fig. 6(g)] is consistently identified as vegetation. 




FIG. 6 (A) "FELTWELL-LHV" SAR IMAGE OF U.K. (309 X 433), (B) 
SEGMENTATION OBTAINED BY THE SWE, (C) SEGMENTATION 
OBTAINED BY KWE, (D) SEGMENTATION OBTAINED BY THE 
PROPOSED KWE-GA, (E) [RESPECTIVELY, (F) AND (G)] ZOOM 
OF AN AREA EXTRACTED FROM (B) [RESPECTIVELY, (C) AND 

(D)] 
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Conclusion 

We have developed a new level set algorithm based on 
kurtosis value of wavelet coefficients energy as a 
feature extractor tools and genetic algorithm as a 
classifier tools for recognition of type of land cover of 
each region of satellite SAR images. A new energy 
named kurtosis wavelet energy (KWE) based on 
genetic algorithm (GA) is proposed to use as a 
recognition tools for texture classification of each 
region in the satellite SAR images. In comparison with 
the KWE single based segmentation model, the KWE 
based on GA achieves better performance on the SAR 
images. It also performs better than the method 
proposed in [15] in often cases. Experimental results 
show that the proposed KWE-GA algorithm is more 
robust for accurate recognition of several land cover 
kinds of satellite SAR images. 
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